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A numerical and analytic treatment is presented here of the evolution of initial data of the kind 
that was conjectured by Hertog, Horowitz and Maeda to lead to a violation of cosmic censorship. 
That initial data is essentially a thick domain wall connecting two regions of anti de Sitter space. 
The evolution results in no violation of cosmic censorship, but rather the formation of a small black 
hole. 
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i I. INTRODUCTION 

o ■ 

A longstanding issue in general relativity is that of cosmic censorship: the question of whether the singularities 
formed in gravitational collapse are hidden inside black holes. This old question was given a new twist by Hertog, 
Horowitz and Maeda 1] who conjectured that cosmic censorship can be violated in certain spacetimes that are asymp- 
totically anti do Sitter. The proposed counterexample consists of a matter model, an initial data set and an argument 
that the evolution of that initial data gives rise to a singularity that cannot be hidden inside a black hole. 
00 The matter model is a scalar field with a potential with two minima below zero. The potential is chosen so that the 

^-H ^ system just barely satisfies the positive mass theorem. The initial data is essentially a spherically symmetric, thick 
domain wall that interpolates between regions of the two anti de Sitter spaces corresponding to the two potential 
minima. The data is chosen so that for an inner region of radius R, the mass is proportional to R. Since anti de Sitter 
space is unstable, it is argued in 1] that the evolution of this initial data will result in a singularity in the central 
region and further that for sufficiently large R the mass of the spacetime will not allow the formation of a black hole 
', large enough to cover the singularity. 
OO ' The nature of this issue was changed when DafermosiJ proved that any singularities formed in the evolution of the 
] system of Q cannot be visible to observers at infinity. Thus if the arguments of were correct, then the singularity 
\l would extend to infinity. In other words the entire space would collapse in a Big Crunch. 

However, various authors UJJ,Jp,j6| have expressed misgivings about the arguments of (ij (including the authors 
O . of in 0)- The authors of [JUQl express doubts over the assumption in that the central region can be well 
O^' approximated as a spacetime that is homogeneous but not anti de Sitter. The authors of |^, |^ raise the possibility 
, that the wall might move outward indefinitely. 

It is not clear how to settle this issue using analytical means. Therefore it makes sense to perform numerical 
^ ' simulations of the evolution of the initial data of Ij and find the outcome. Such a simulation was reported in 3] for 
■ R — 7 with the result that a small black hole formed rather than a naked singularity. However, this did not resolve 
/\ ' the issue since it was argued in^jQ that "sufficiently large" R means R > 600. 

In this work we simulate the evolution of the initial data ofQ for large R. Here too, the result is the formation of 
a small black hole, not a naked singularity. Though the results are numerical, the most important properties of this 
system can be understood analytically using the properties of perturbations of anti de Sitter space. 

Sec. II presents the relevant equations and numerical method. Sec. Ill contains a treatment of perturbations of 
anti de Sitter space. Results of the numerical simulations are presented in Sec. IV and conclusions in Sec. V. 



II. EQUATIONS AND NUMERICAL METHODS 



The system to be studied is a spherically symmetric scalar field (f> with a potential V. The appropriate equations 
are therefore the Einstein-scalar equations: 

Gab = V,0Vfa0-ffah(iV^0Ve(/.+ l/) (1) 

dV 

^-^"^ = IJT (2) 
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(Here we are using units where SttG = 1). We use polar-radial coordinates for the metric which puts it in the form 

ds^ = -a^dt^ + a^dr^ + {dO^ + sij^OdLp^) (3) 

It is helpful to define the quantities X and Y given hy X = d4)/dr and Y = {a/ a)d(f)/ dt. Then equation Q yields 
the following evolution equation for Y 

dY I d ( dV 

a^^^a^i- (4) 

while Einstein's equation yields the following equations which are used to find the metric components a and a. 

da ail — i / o o 9 % 

- ^ ' ' \ra{X^ + Y^ + 2a^V) (5) 



dr 2r 



l-ln(aa) = '-(X^ + Y^) (6) 

Einstein's equation also provides an evolution equation for a. Defining the quantity C = da/dt — raXY/2 then it 
follows from Einstein's equation that C vanishes. We will use this as a code check. That is, equations H4I6|) are used 
to simulate the evolution of the system and then the results of that simulation are checked by seeing that C converges 
to zero. 

The numerical method used is essentially that of reference Q. That is spatial derivatives are replaced by cen- 
tered differences for unequally spaced points, while time evolution is done using the iterated Crank-Nicholson (ICN) 
method. "s] The equations are stabilized using Kreiss-Oliger dissipation.'^ 

The spacing of the grid points is different from that of 3] . Define the coordinate p by r = tan p. Though the 
coordinate r is used in all equations, we choose a grid of points that is equally spaced in p. 

The potential used is that of reference^ 

V{(l)) = -3 + 50(/>2 - 8103 f^^6 

where the constant k is chosen so that the system just barely satisfies the positive mass theorem. The initial data 
is chosen to minimize the contribution of the potential to the total mass. Here the minimum is found over all field 
configurations that are in the true vacuum at r = and the false vacuum for r > R. This leads to an ordinary 
differential equation for (j) that is solved using a shooting method as described inQ- For our purposes, an important 
property of the solution is that in the central region (/> oc where /3 = (-\/409 — 3)/2 « 8.6. 



III. ANTI DE SITTER PERTURBATIONS 



We now consider the behavior of the central region as a perturbation of anti de Sitter space. The initial data in the 
central region has (f> oc where f3 ~ 8.6. Thus the scalar field is very small. Keeping only terms up to linear order 
in (p we find that equations (|5I6|I yield — a^^ — 1 + r^, while equation 10} becomes 



r2[l + r2]-^ -100(1 + r2)0 = O (8) 



Defining tp = np and p = tan~"'^r we obtain 

dt"^ dp^ cos^/9 

Provided that the spatial and temporal variation of ip is sufficiently large, the last term in this equation can be 
neglected yielding for an ingoing wave ^ f{t + p) or 

(/) = i/(t + tan-V) (10) 
r 

We will later see from the numerical results that for large walls, the last term in equation @ is negligible. However, 
the reason for this is easy to understand analytically: Define a "coordinate thickness" Ar of the wall as follows: Let 
Ar be the amount that r varies as 4> varies from 10% to 90% of its maximum value. Correspondingly define Ap. Since 
the initial data depends only on r/R, it follows that Ar cx R. However, this means that for large walls {i.e. those 
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FIG. 1: the scalar field </> at i = 0. i? = 100 and = 256000 
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FIG. 2: the constraint at t = 0.5 for R = 100 at two different resolutions: C for iV = 256000 (sohd line) and 8C for N = 512000 

(dashed line) 

with i? » 1) we have Ap oc R~^. Thus viewed as a function of p, the initial data is essentially zero in the central 
region and then steeply increases to its maximum value in a narrow region near 7r/2. The evolution of such initial 
data is a narrow wave packet that propagates inward. 

Once the wave packet reaches values of r <i; 1, the difference between anti de Sitter space and Minkowski space 
becomes irrelevant. The subsequent behavior of the wave packet then depends on whether it will shrink to a size 
smaller than its Schwarzschild radius or whether it will disperse. 

IV. NUMERICAL RESULTS 

All runs were done in double precision on a Sun Blade 2000. Let iV + 1 be the number of spatial grid points. We 
first present results for R = 100. Figure 1 shows the initial data for The evolution of this initial data produces 
very narrow wave packets; so all plots of the evolution of the initial data will be confined to the range of r where the 
wave packet is present. 

Figure 2 shows the result of a convergence test of the code. Here, the system is evolved to a time of 0.5 in two 
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FIG. 3: the scalar field cj) at t = 0.731. R = 100 and N = 256000 
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FIG. 4: y at t = 0.731. R = 100 and iV = 256000 

different runs: one with N = 256000 and the other with N — 512000. What is plotted in the figure is the quantity C 
for the N = 256000 run (solid line) and 8C for the N = 512000 run (dashed line). The result shows that the constraint 
converges to zero at third order. This is a bit surprising, since from the finite difference techniques used one would 
expect the code to be only second order convergent. However, it may be that for this type of ingoing wavepacket the 
leading part of the truncation error for the constraint vanishes. 

Figures 3-9 show the results of a run with TV ~ 256000. Figures 3 and 4 show respectively (j) and Y for the evolution 
of the initial data to a time t=0.731. Note that d(t)/dt ^ 5700. Therefore, the last term in equation is negligible, 
justifying the expression in equation (|10|) . Figure 5 shows r0 as a function oit + p for two times: t ~ 0.2 (solid line) 
and t = 0.4 (dashed line). These curves agree, showing that equation (|10|1 is a good approximation for this part of 
the evolution. In contrast, figure 6 shows (/) as a function of r for these two times. Here one sees two well separated 
wave packets. 

Figures 7, 8 and 9 show respectively (/), Y and a at the final time. Here the final time, which is slightly less than 
7r/2, is chosen by having the simulation end when the maximum value of a reaches 3.5 corresponding to 2m/r — 0.92. 
The reason for this is that the coordinate system used here breaks down when a trapped surface forms and that 
breakdown is signalled by a becoming large where the trapped surface forms. Thus figure 9 indicates that a trapped 
surface forms at r « 0.0056, while figures 7 and 8 indicate that the region of high curvature is contained within this 
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FIG. 5: r(j> as a function oi t + p at t = 0.2 (solid line) and t = 0.4 (dashed line). R = 100 and N = 256000 
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FIG. 6: as a function of r at f = 0.2 (solid line) and t = 0.4 (dashed line). R = 100 and iV = 256000 



trapped surface. Thus the resuh of the evolution of the initial data with R = 100 is a small black hole. 

These numerical results for R = 100 do not directly address the conjecture of [l| which is for R > 600. However, 
a wall with R ~ 600 leads to a wavepacket that is much more narrow in p and therefore requires far more resolution 
and thus much more computer memory and time for a simulation that evolves all the way to black hole formation. 
Instead we will evolve such initial data for a comparatively short time that is nonetheless long enough to (i) verify 
that the result is a narrow wavepacket and (ii) show that the wavepacket has enough energy that its collapse will 
result in a black hole. Figures 10 and 11 show the results of a run with R = 600 and N = 1024000 evolved to a time 
of t = 0.1. Here (p is plotted in figure 10 and Y is plotted in figure 11. Note that the result is a narrow wavepacket 
centered near r = 10, just as one would expect from the treatment of section III. Also note that d<j)/dt ^ 3.8 x 10^ (j). 
For comparison, figures 12 and 13 give respectively and Y for a run with R — 100 and t = 0.1. {N — 256000 
for this run). Note that the R — 600 wavepacket has a higher amplitude and shorter wavelength than the i? = 100 
wavepacket. Thus the R = 600 wavepacket has more energy within a shell of smaller thickness than the i? = 100 
wavepacket. Since a black hole forms from the further evolution of the R = 100 case, it then follows that the evolution 
of the more energetic wavepacket of the R = 600 case will also result in the formation of a black hole. 
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FIG. 7: at the final time. R = 100 and iV = 256000 
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FIG. 8: Y at the final time. R = 100 and iV = 256000 



V. CONCLUSIONS 



We then see that the initial data of [ij when evolved, does not form a naked singularity, but instead forms a 
small black hole. In hindsight, the reason for this is clear. The authors of [j| assumed, correctly, that the evolution 
of the central region could be described for a while by a perturbation of anti de Sitter space. However, they also 
assumed, incorrectly, that the perturbation is homogeneous. In the initial data (j) oc where (3 ~ 8.6. Thus the initial 
amplitude for the homogeneous part of the perturbation is so small as to be negligible. Instead, the perturbation 
comes from the wall, i.e. the transition region between the two vacua. In the p coordinate, this transition region is 
very narrow (especially for walls with large R) . Therefore the evolution of this initial data is a narrow wavepacket 
that propagates inward. In the course of the evolution, the wavepacket becomes sufficiently concentrated that it is no 
longer well described by perturbation theory. The numerical simulation then shows that the wavepacket collapses to 
form a small black hole. Thus cosmic censorship is not violated by the evolution of the initial data given in 
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FIG. 10: (j>a.tt = 0.1. R = 600 and N = 1024000 
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